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Abstract 

Background: Tularemia is a zoonotic disease caused by Francisella tularensis that has been found in many different 
vertebrates. In Germany most human infections are caused by contact with infected European brown hares {Lepus 
europaeus). The aim of this study was to elucidate the epidemiology of tularemia in hares using phenotypic and 
genotypic characteristics of F tularensis. 

Results: Cultivation of F tularensis subsp. holarctica bacteria from organ material was successful in 31 of 52 hares 
that had a positive PCR result targeting the Ft-M19 locus. 17 isolates were sensitive to erythromycin and 14 were 
resistant. Analysis of VNTR loci (Ft-M3, Ft-M6 and Ft-M24), INDELs (Ftind33, Ftind38, Ftind49, RD23) and SNPs 
(B.17, B.I 8, B.19, and B.20) was shown to be useful to investigate the genetic relatedness of Francisella strains in this 
set of strains. The 14 erythromycin resistant isolates were assigned to clade B.I, and 16 erythromycin sensitive 
isolates to clade B.IV and one isolate was found to belong to clade B.ll. MALDI-TOF mass spectrometry (MS) was 
useful to discriminate strains to the subspecies level. 

Conclusions: F. tularensis seems to be a re-emerging pathogen in Germany. The pathogen can easily be identified 
using PCR assays. Isolates can also be identified within one hour using MALDI-TOF MS in laboratories where specific 
PCR assays are not established. Further analysis of strains requires genotyping tools. The results from this study 
indicate a geographical segregation of the phylogenetic clade B.I and B.IV, where B.I strains localize primarily within 
eastern Germany and B.IV strains within western Germany. This phylogeographical pattern coincides with the 
distribution of biovar I (erythromycin sensitive) and biovar II (erythromycin resistance) strains. When time and costs 
are limiting parameters small numbers of isolates can be analysed using PCR assays combined with DNA 
sequencing with a focus on genetic loci that are most likely discriminatory among strains found in a specific area. 
In perspective, whole genome data will have to be investigated especially when terrorist attack strains need to be 
tracked to their genetic and geographical sources. 
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Background 

Tularemia is a rare zoonotic disease caused by Francisella 
tularensis, a Gram negative, facultative intracellular, fas- 
tidious bacterium [1]. Most infections in animals and 
humans are caused by two F. tularensis subspecies, F. 
tularensis subsp. tularensis (Jellison type A) and F. 
tularensis subsp. holarctica (Jellison type B). F. tularensis 
type A is endemic in North America and type B is located 
in Europe, Asia, and North America [2-4]. Three biotypes 
of the less virulent type B have been described: biovar I 
(erythromycin sensitive), biovar II (erythromycin resist- 
ant), and biovar japonica which can ferment glycerol [4]. 

In Germany, human infections are usually caused by 
skinning, preparing or eating infected hares or drinking 
contaminated water. F. tularensis was sporadically diag- 
nosed in humans in the first half of the 20th century in 
Germany but almost disappeared in the following decades 
[5,6]. Between 1983 and 1992 only four sporadic cases of 
tularemia were notified in hares or rabbits from Lower 
Saxony, Rhineland-Palatinate, North Rhine-Westphalia and 
Baden- Wurttemberg, respectively [6]. After years without 
reported cases in animals the re-emergence of tularemia 
started in 2004 with an outbreak of tularemia in a semi-free 
living group of marmosets (Callithrix jacchus) in Lower 
Saxony [7], and in December 2005 an outbreak with 15 hu- 
man cases due to contact with infected hares was reported 
from Hesse [8]. The detection of F. tularensis subsp. 
holarctica in organ samples of these hares using PCR assays 
was the beginning of our investigations of tularemia in 
European brown hares (Lepus europaeus) in Germany. 

A variety of PCR methods has been established for the 
detection of F. tularensis DNA in both clinical and envir- 
onmental specimens [9-11]. Farlow et al. developed a typ- 
ing assay based on the variable-number of tandem repeats 
(VNTRs) [12] and Johansson et al. also described a 
twenty-five VNTR marker typing system that was used to 
determine the worldwide genetic relationship among F. 
tularensis isolates [1]. Bystrom et al. selected six of these 
25 markers that were highly discriminatory in a study of 
tularemia in Denmark [13]. Vogler et al. [14] investigated 
the phylogeography of F. tularensis in an extensive study 
based on whole-genome single nucleotide polymorphism 
(SNP) analysis. From almost 30,000 SNPs identified 
among 13 whole genomes 23 clade- and subclade-specific 
canonical SNPs were identified and used to genotype 496 
isolates. This study was expanded upon in another study 
that used a combination of insertion/deletions (INDELs) 
and single nucleotide polymorphism analysis [15]. 

The aim of this study was to elucidate the molecular epi- 
demiology of F. tularensis in European brown hares in 
Germany between 2005 and 2010. Several previously pub- 
lished typing markers were selected and combined in a 
pragmatic approach to test whether they are suitable to elu- 
cidate the spread of tularemia in Germany. This included 



cultivation, susceptibility testing to erythromycin, a PCR 
assay for subspecies differentiation detecting a 30 bp dele- 
tion in the Ft-M19 locus, VNTR typing, INDEL, SNP, and 
MALDI-TOF analysis. This is important because it im- 
proves our understanding of the spread of tularemia and 
may help to recognize outbreaks that are not of natural 
origin. 

Results 

Cultivation and identification of isolates 

Cultivation of bacteria from organ specimens was suc- 
cessful in 31 of 52 hares which had a positive PCR result 
targeting the locus Ft-M19 that was also used to differ- 
entiate F. tularensis subsp. holarctica from other F. 
tularensis subsp. [11]. F. tularensis subsp. holarctica was 
identified in all 52 cases. 

Biovars 

Seventeen isolates were susceptible to erythromycin cor- 
responding to biovar I, whereas fourteen were resistant 
(biovar II). The geographic distribution is given in 
Table 1, Figure 1 and the susceptibility of the isolates in 
Additional file 1: Table S2. 

VNTR typing 

In a pilot study, six loci (Ft-M3, Ft-M6, Ft-M20, Ft-M21, 
Ft-M22, and Ft-M24) were amplified and sequenced, but 
only the loci Ft-M3, Ft-M6, and Ft-M24 were discrimin- 
atory. The strains tested in the pilot study are indicated 
in Additional file 1: Table S2(*). The following identical 
results were obtained for all these strains: Ft-M20: 255 
bp; Ft-M21: 403 bp; Ft-M22: 241 bp. The loci Ft-M3 and 
Ft-M6 (repeat: TTG GTG AAC TTT CTT GCT CTT) 
were further used to analyse DNA samples extracted 
from cultivated bacteria. Sequencing of Ft-M3 identified 
two different repetitive elements, Ft-M3a (ATC CTT 
ATT), and Ft-M3b (GTC TTT GTT), respectively. The 
number of these repeats was determined separately. The 
size obtained for Ft-M24 (repeat: ATA AAT TAT TTA 
TTT TGA TTA) correlated with the size observed previ- 
ously for the B.IV (B.18) clade. All VNTR results are 
given in Additional file 1: Table S2. 

INDEL analysis 

Conventional PCR assays with subsequent gel electro- 
phoretic size determination of the amplicon allowed clear 
discrimination between amplicons with or without the re- 
spective deletions which was confirmed by sequencing in 
some cases (data not shown). The 31 isolates showed four 
different INDEL patterns (Additional file 1: Table S2). 
Based on INDELs and SNPs (see below) 14 isolates were 
assigned to clade B.I (B.20), one to B.II (B.17), and 16 to 
B.IV (B.18), according to the nomenclature in Karlsson 
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Table 1 Original and geographic data of Francisella tularensis subsp. holarctica isolates (BW - Baden-Wurttemberg, 
BY - Bavaria, NRW - North Rhine-Westphalia, LS - Lower Saxony, SN - Saxony, TH - Thuringia) 



Year 


Strain number 


Site 


Federal state 


Latitude [°NORTH] 


Longitude [°EAST] 


Altitude [m] 


2006 


06T0001 


Moorgrund 


TH 


50,838005 


10,291767 


279 


2007 


08T0001 


Dingelstadt 


TH 


51,315205 


10319329 


335 


2007 


08T0008 


Allersberg 


BY 


49,251389 


11,234261 


388 


2007 


08T0010 


Sehnde 


LS 


52,31262 


9,967105 


71 


2007 


08T0013 


Ehingen 


BY 


49,300734 


10,571476 


415 


2008 


08T0014 


Weissach-Flacht 


BW 


48,833991 


8,91309 


406 


2008 


08T0015 


Leonberg-Hofingen 


BW 


48,816676 


9,016877 


379 


2007 


08T0070 


Einbeck-Kohnsen 


LS 


51,707717 


10,000538 


121 


2008 


08T0071 


Brake 


LS 


53,326329 


8,478167 


2 


2008 


08T0072 


Gottingen-Roringen 


LS 


51,532638 


9,92816 


153 


2008 


08T0073 


Twulpstedt-Rummer 


LS 


52,224403 


11,01102 


127 


2008 


08T0075 


Wurzburg 


BY 


49,794256 


9,927489 


173 


2009 


09T0105 


Geseke 


NRW 


51,639416 


8,509738 


105 


2009 


09T0108 


Geseke 


NRW 


51,639416 


8,509738 


105 


2009 


09T0109 


Geseke 


NRW 


47,724358 


9,406025 


518 


2009 


09T0114 


Markdorf 


BW 


51,639416 


8,509738 


105 


2009 


09T0116 


Geseke 


NRW 


51,444502 


12,169177 


115 


2009 


09T0146 


Wiedemar 


SN 


48,864962 


9,024819 


337 


2008 


09T0163 


Hemmingen 


BW 


53,770141 


7,693722 


1 


2008 


09T0164 


Spiekeroog 


LS 


53,770141 


7,693722 


1 


2008 


09T0165 


Spiekeroog 


LS 


53,770141 


7,693722 


1 


2008 


09T0166 


Spiekeroog 


LS 


53,770141 


7,693722 


1 


2008 


09T0167 


Spiekeroog 


LS 


53,770141 


7,693722 


1 


2008 


09T0169 


Spiekeroog 


LS 


53,745892 


7,480842 


4 


2008 


09T0171 


Langeoog 


LS 


53,745892 


7,480842 


4 


2009 


09T0179 


Langeoog 


LS 


53,745892 


7,480842 


4 


2010 


10T0014 


Hemmingen 


BW 


52,314054 


9,722783 


56 


2010 


10T0115 


Waltrop 


NRW 


51,624087 


7,39465 


69 


2010 


10T0125 


Geseke 


NRW 


51,639223 


8,469223 


103 


2010 


10T0128 


Empfingen 


BW 


48,391957 


8,708282 


491 


2010 


10T0131 


Oppenweiler 


BW 


48,985678 


9,460219 


270 



et al 2013 [16], where the B.I clade was re-defined to in- 
clude both Bl and B3 in Svensson et al [15]. 



SNP typing 

The results of SNP typing are given in Additional 
file 1: Table S2. All strains used in a pilot study 
showed ancestral states of the following SNPs: B.16: 
G; B.21:G; B.22:G; B.24:T. These strains are indicated 
(*) in Additional file 1: Table S2. Therefore, only the 
SNPs B.17, B.18, B.19, and B.20 were further investi- 
gated for all isolates. 



MALDI-TOF MS analysis 

All isolates (n=31) yielded high quality spectra. MALDI- 
TOF was found to be useful for rapid identification of 
isolates to subspecies level within one hour. However, 
the obtained clusters (Figure 2) did not conform to the 
genetic clusters (Additional file 1: Table S2). 

Geographical clustering 

Cases of tularemia in hares were identified in eight of six- 
teen federal states of Germany reaching from islands in the 
North Sea to regions at Lake Constance in the southern 
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Figure 1 Germany: areas where Francisella tularensis subsp. holarctica isolates were found in hares. Erythromycin sensitive strains occur in 
regions marked with blue dots, erythromycin resistant strains occur in regions marked with red dots. 



part of Germany. All cases were found below 500 m above 
sea level. Isolates belonging to biovar I could be found in 
the western part of Germany whereas biovar II occurred in 
the eastern region (Table 1 and Additional file 1: Table S2, 
Figure 1). Molecular typing resulted in further discrimin- 
ation of clusters within the biovars. Isolates resistant to 
erythromycin and genetically assigned to clade B.I were 
found only in Lower Saxony, Thuringia, Bavaria and Sax- 
ony. Strains that were sensitive to erythromycin could be 
assigned to clade B.II (Ftind38) and B.IV (B.18) as given in 
Additional file 1: Table S2. 

Stability testing 

The investigated markers for two Francisella isolates 
(06T0001 from hare and 10T0191 from fox) were stable 
even after 20 passages in cell culture and had identical 
results for the markers Ft-M3 (297 bp), Ft-M6 (311 bp), 



Ftind33 (deletion), Ftind38 (insertion), and Ftind49 
(insertion). 

Discussion 

In Thuringia the first case of tularemia in a hare was 
reported in 2006 [17]. In Lower Saxony 2,162 European 
brown hares and European rabbits (Oryctolagus cuniculus) 
were screened for tularemia between 2006 and 2009 using 
cultivation and PCR assays. Francisella specific PCR as- 
says were positive in 23 hares and 1 rabbit which were fur- 
ther confirmed by cultivation of F. tularensis subsp. 
holarctica in 12 hares [18]. In the present study, cases of 
tularemia in hares in Germany from 2005 to 2010 were in- 
vestigated. During this period a total of 52 hares were 
found positive in PCR assays for F. tularensis subsp. 
holarctica DNA and from 31 of these cases Francisella 
strains could be isolated. MALDI-TOF analysis was also 
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Figure 2 Dendrogram constructed from MALDI-TOF mass spectrometry spectra of 31 Francisella tularensis ssp. holarctica strains and 
representatives of ssp. tularensis, mediasiatica, and novicida. 



used to rapidly identify Francisella to the subspecies level 
as was previously shown by Seibold et al [19]. 

Several positive specimens were found on the North Sea 
islands Langeoog and Spiekeroog (LS), around Soest (NR), 
Darmstadt (H), and Boblingen (BW). These natural foci 
and also sporadic cases in other regions of Germany were 
found below 500 m above sea level In the Czech Republic 
typical natural foci of tularemia occurred in alluvial forests 
and field biotopes below 200 m sea level with mean an- 
nual air temperature between 8.1-10.0°C and mean annual 
precipitation of 450-700 mm [20]. In Germany, an out- 
break of tularemia in a colony of semi-free living marmo- 
sets was located in a region with geographic and 
ecological conditions similar to the hare habitats in the 
Czech Republic: field biotopes 175 m above sea level 
(<200 m) with 9.2°C mean annual air temperature and 
642 mm mean annual precipitation [8]. In Germany, tular- 
emia of hares occurs in regions with rather humid soil like 
in alluvial forests and alongside rivers, but this obviously 
corresponds with the natural habitat of hares. 

Specimens were screened using a PCR assay targeting 
Ft-M19 described by Johansson et al. [11] which allows 
the simultaneous identification of the species F. tularensis 
and the differentiation of the subspecies holarctica from 
other (sub-) species. All samples could be attributed to 
F. tularensis subsp. holarctica. 

We found a clear segregation of clade B.I and clade B.IV 
in Germany, B.I strains dominate in eastern Germany and 
B.IV within western Germany (Figure 1). Clade B.I is 
known to dominate in Europe between Scandinavia and 
the Black Sea [15,16,21-23]. The other dominating 



European clade is B.IV (B.18) which can be found over a 
large area of western and central Europe, and, based upon 
this study, western Germany [21,23-26]. We found only 
one strain of the B.II clade isolated in Bavaria. Strains of 
the B.II clade are most frequently isolated in the USA, but 
are found sporadically in Europe as well [16,21]. 

The phylogeographical pattern of clade B.I and B.IV, co- 
incide with the geographical distribution of biovar II and 
biovar I strains, respectively. Previously, biovar I strains 
(erythromycin sensitive) have been reported from Western 
Europe (France, Germany, Spain and Switzerland), North- 
America, Eastern Siberia and the Far East while biovar II is 
present in the European part of Russia as well as Northern, 
Central and Eastern Europe (Austria, Germany, Sweden 
and Turkey) [27-31]. A mixture of both biotypes has been 
reported in Sweden, Norway, Bulgaria, Russia and 
Kazakhstan [27,28,32]. Isolation of both biovars from ro- 
dents in a single settlement in Moscow as well as from 
water samples collected in the Novgorod region [27] indi- 
cate coexistence of the biovars in the same epidemio- 
logical foci. Taken together, a geographical separation of 
F. tularensis strains seems to exist in Germany. The 
phenotypically defined biovar I (erythromycin sensitive) 
and phylogenetically defined clade B.IV strains are con- 
fined in western Germany, whereas biovar II (erythro- 
mycin resistance) and clade B.I strains cluster in eastern 
Germany. This is interesting and may reflect a competi- 
tion between the two subpopulations or unknown 
underlying ecological or epidemiological differences. 

A deletion in the genome of F. tularensis subsp. 
holarctica in RD 2 3 is typical for strains of F. tularensis 
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subsp. holarctica in France, the Iberian Peninsula and also 
Switzerland, where biovar I predominates [24,25,27]. 
However, in one erythromycin susceptible isolate from 
Bavaria (08T0013), classified in this study as belonging to 
the B.II clade, RD 2 3 was not deleted, thus showing that de- 
letion of RD 2 3 is not correlated with sensitivity to erythro- 
mycin. The molecular mechanisms of resistance to 
erythromycin have not been functionally established, but 
mutations identified in domain V of the 23S rRNA of 
biovar II strains, could provide a likely explanation [33]. 

Although 25 VNTR markers have been described for 
the typing of Francisella, it is pragmatic to investigate only 
loci of interest depending on the prevalent subspecies of 
F. tularensis, the efficiency of PCR assays for single loci, 
and existing data [1,13,34]. Sequence analysis of the locus 
Ft-M3 resulted in two different repeats denominated here 
as Ft-M3a corresponding with SSTR9E and Ft-M3b corre- 
sponding with SSTR9A as described previously by 
Johansson et al. [35]. Johansson et al. and Bystrom et al. 
also found that locus Ft-M3 is the most variable marker 
[1,13]. In the Francisella genome variations of DNA se- 
quences in spite of identical repeat length have been de- 
scribed for short-sequence tandem repeats [35,36]. Locus 
Ft-M6 showed less variability with only three PCR frag- 
ment sizes being observed among the strains. We obtained 
the same amplicon sizes that were described in previous 
studies for locus Ft-M3 (Additional file 1: Table S2) 
[14,37] and for locus Ft-M6 (Additional file 1: Table S2) 
[14,37]. Svensson et al. developed a sophisticated real-time 
PCR array for hierarchical identification of Francisella iso- 
lates [15]. Only three (Ftind33, Ftind38, Ftind49) out of 
five INDEL loci were discriminatory among our set of 
F. tularensis subsp. holarctica isolates. Ftind48 is a marker 
for B.I to B.IV clades (non-japonica/non-california) and is 
not expected to vary for these isolates, and Ftind50 is 
targeting a specific deletion that so far only has been 
found in LVS. It was possible to simplify these assays to 
conventional PCR assays that allowed a simple read out 
based on gel electrophoresis. We identified clusters of 
strains that had the same INDELs and SNPs as strains de- 
scribed by Svensson et al. [15]. In our study the analysis of 
VNTR and INDELs of two F. tularensis subsp. holarctica 
strains (06T0001, 10T0191) that were passaged twenty 
times in Ma- 104 cells showed that these genomic ele- 
ments were stable. Johansson et al. demonstrated for two 
VNTR loci (SSTR9 and SSTR16) that they were actually 
stable over 55 passages [35]. The VNTR pattern for strains 
belonging to clade B.I was more variable compared with 
the pattern obtained for clade B. IV (Additional file 1: 
Table S2), as was observed previously [21,23-25]. This 
might indicate that clade B.IV is more recently introduced 
in Germany than clade B.I. 

We have applied several typing tools in a polyphasic ap- 
proach in order to determine their value for identifying 



groups of Francisella strains in Germany. We found 
strains belonging to biovars I and II of F. tularensis subsp. 
holarctica. Although SNP loci are the most informative 
markers for typing of Francisella this method may have to 
be adapted to local strains [37,38]. 

Conclusions 

F. tularensis seems to be a re-emerging pathogen in 
Germany that infects hares in many regions and causes 
a potential risk for exposed humans such as hunters and 
others who process hares. The pathogen can easily be 
identified using PCR assays directly on DNA extracted 
from organ specimens or cultivated strains. Isolates can 
also be identified rapidly using MALDI-TOF MS in 
routine laboratories where specific PCR assays for 
F. tularensis are not established. To identify differences 
and genetic relatedness of Francisella strains, analysis of 
VNTR loci (Ft-M3, Ft-M6 and Ft-M24), INDELs 
(Ftind33, Ftind38, Ftind49, RD23) and SNPs (B.17, B.18, 
B.19, and B.20) was shown to be useful in this set of 
strains. When time and costs are limiting parameters 
isolates can be analysed using simplified PCR assays with 
a focus on genetic loci that are most likely discrimin- 
atory among strains found in a specific area. For the fu- 
ture whole genome sequencing using next generation 
sequencing is desirable and should provide more genetic 
information of Francisella strains. Based on these data a 
more detailed view on the epidemiology of tularemia will 
become possible [39]. 

Methods 

Samples 

Organ specimens (e.g. spleen, liver, lung, and/or kidney) 
of European brown hares that were suspicious of tular- 
emia were collected by local veterinary authorities in 
Germany since 2005 and sent for confirmatory testing to 
the National Reference Laboratory for Tularemia of the 
Friedrich-LoefQer-Institut in Jena. Francisella strains 
were cultivated on cysteine heart agar (Becton Dickinson 
GmbH, Heidelberg, Germany) supplemented with 10% 
chocolatized sheep blood and antibiotics in order to sup- 
press the growth of contaminants. One litre of culture 
medium contained 100 mg ampicillin ( Sigma- Aldrich 
Chemie, Taufkirchen, Germany) and 600 000 U poly- 
myxin B (Sigma- Aldrich Chemie). Plates were incubated 
at 37°C with 5% C0 2 for up to 10 days. Typical colonies 
are grey-green, mostly confluent, glossy, and opaque. 
Gram staining was performed routinely and showed 
Gram negative coccoid bacteria. 

The reference strains F. tularensis subsp. tularensis 
(FSC 237), mediasiatica (FSC 147), and F. novicida 
(ATCC 15482) were obtained from the Bundeswehr 
Institute of Microbiology, Munich, Germany, and F. 
philomiragia (DSMZ 7535) was obtained from the 
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German Collection of Microorganisms and Cell Cul- 
tures, Braunschweig, Germany, respectively. 

Erythromycin susceptibility 

All F. tularensis subsp. holarctica isolates were tested for 
their erythromycin susceptibility using Erythromycin discs 
[30 ug] and M.I.C.Evaluator™ (Oxoid, Wesel, Germany) in 
order to discriminate the susceptible biovar I from the re- 
sistant biovar II as described previously [40]. 

DNA extraction 

50 mg of organ material or 200 ul of cell culture super- 
natant were lysed and the DNA was extracted using the 
High Pure PCR Template Preparation Kit (Roche Diag- 
nostics, Mannheim, Germany) according to the manufac- 
turer's instructions. If cultivation was successful some 
colonies were resuspended in 200 ul phosphate-buffered 
saline, boiled at 90°C for 10 minutes and DNA was pre- 
pared as described above. Finally, DNA was eluted in 200 
ul elution buffer. 5 ul were applied in each PCR assay. 

Diagnostic PCR assay 

F. tularensis subsp. holarctica was identified using a 
PCR assay with primer pair C1/C4 targeting the locus 
Ft-M19 that distinguishes the two major subspecies F. 
tularensis subsp. holarctica and F. tularensis subsp. 
tularensis which was carried out as described by 
Johansson et al. [11]. 

VNTR typing 

In pilot experiments 6 VNTR loci (Ft-M3, Ft-M6, Ft-M20, 
Ft-M21, Ft-M22, and Ft-M24) were investigated as de- 
scribed by Bystrom et al. [13]. The loci found discrimin- 
atory were then subsequently analysed in all 31 isolates. 
The amplification of the VNTR loci was carried out under 
the same cycling conditions as the diagnostic PCR assay 
except for the annealing temperature of 56°C. The frag- 
ments were cut out of the agarose gel and DNA was puri- 
fied using the innuPrep Gel Extraction Kit (Analytik Jena 
AG, Jena, Germany) according to the manufacturers 
instructions. Subsequently, DNA amplificates were se- 
quenced as described below. 

INDEL analysis 

Five INDELs (Ftind33, Ftind38, Ftind48, Ftind49, and 
Ftind50) that are discriminatory among F. tularensis 
subsp. holarctica were selected from the loci described 
by Svensson et al. [15]. The real-time PCR assays with 
melting curve analyses were simplified by using conven- 
tional PCR assays. The primers "CP" and "OUT" for the 
respective loci were used as described by Svensson et al. 
The reaction mixture consisted of 5 ul 10 x PCR buffer 
with 1.5 mM MgCl 2 (Genaxxon, Stafflangen, Germany), 
2 \A of dNTP mix (each 2 mM, Carl Roth GmbH, 



Karlsruhe, Germany), 1 \A of each primer, 0.2 \A of Taq 
DNA polymerase (5 U/ul, Genaxxon), 5 ul of DNA ex- 
tract and deionised water to a final volume of 50 ul. 
After denaturation at 95°C for 5 min, 35 cycles of ampli- 
fication were performed with denaturation at 95°C for 
30 s, primer annealing at 60°C for 60 s, and primer ex- 
tension at 72°C for 30 s. After a final extension step at 
72°C for 30 s amplicons were separated using 2.5% agar- 
ose gel electrophoresis and visualized using ethidium 
bromide staining under UV light. 

SNP typing 

Four often SNPs (B.17, B.18, B.19, and B.20) that have been 
found to be useful for the typing of F. tularensis subsp. 
holarctica strains were selected from the loci described by 
Svensson et al. [15]. The primers "C" and "D" for the 
respective loci described by Svensson et al. were used, but 
the primers "D" were shortened by removing the SNP spe- 
cific last nucleotide and the non-binding GC-rich tails that 
were originally added to the allele-specific primer (i.e. 
gcgggcagggcggc). SNPs were detected by sequence analysis 
of the PCR products. The nomenclature used for the desig- 
nation of clades is according to Karlsson et al. 2013 [16]. 

DNA sequencing 

Purified DNA fragments were subjected to cycle sequen- 
cing with BigDye™ Terminator Cycle Sequencing Ready 
Reaction Kit (Applied Biosystems, Darmstadt, Germany). 
Amplification primers were also used as sequencing 
primers. Nucleotide sequences were determined on an 
ABI Prism 310 Genetic Analyzer (Applied Biosystems). 

Analysis of sequence data 

VNTR sequence data were aligned using BioEdit 
(Biological sequence alignment editor, Ibis Therapeutics, 
Carlsbad, CA, USA). 

Stability testing 

The stability of the markers Ft-M3, Ft-M6, Ftind33, 
Ftind38, and Ftind49 was assessed for two F. tularensis 
subsp. holarctica strains that were isolated from a hare 
(06T0001) and a red fox (Vulpes vulpes) (10T0191), re- 
spectively. The isolates were passaged twenty times on 
MA-104 cells in 12.5 ml cell culture flasks (Becton 
Dickinson GmbH, Heidelberg, Germany). Confluent mono- 
layers of MA-104 cells were washed with phosphate- buff- 
ered saline, pH 7.4. The bacterial suspensions or cell 
culture samples were inoculated on the cells at 37°C for 
1 h. The inoculum was replaced with Dulbecco s Modified 
Eagles Medium (DMEM) and incubated at 37°C in a hu- 
midified air atmosphere with 5% C0 2 . After incubation for 
3 to 5 days when the cells detached from the surface, the 
bacteria were harvested by two freeze-thaw cycles. The bac- 
teria/cell suspensions were used for preparation of DNA. 
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MALDI-TOF typing 

Samples were taken from single colonies, ethanol- 
precipitated and extracted with 70% formic acid as de- 
scribed by Sauer et al [41]. The extract was diluted with 
one volume acetonitrile and 1.5 \iL of the mixture was 
spotted to a steel MALDI target. The dried extract was 
overlaid with 1.5 \iL of a saturated solution of 
a-cyano-4-hydroxycinnamic acid in 50% acetonitrile/ 
2.5% trifluoroacetic acid as matrix and was again allowed 
to dry. A custom-made database of reference spectra, or 
main spectra (MSP), was constructed using the BioTyper 
software (version 1.1, Bruker Daltonics, Bremen, 
Germany) following the guidelines of the manufacturer. 
Each sample was spotted six-fold and four single spectra 
with 500 laser pulses each were acquired from each spot 
with an Ultraflex I instrument (Bruker Daltonics) in the 
linear positive mode in the range of 2,000 to 15,000 Da. 
Acceleration voltage was 25 kV and the instrument was 
calibrated in the range of 4,364 to 10,299 Da with refer- 
ence masses of an extract of an Escherichia coli DH5- 
.alpha; strain prepared according to Sauer et al. [41]. 
MSP were generated within the mass range of 2,500 to 
15,000 Da with the following default parameters: com- 
pression of the spectrum data by a factor of 10, baseline 
smoothing by the Savitsky-Golay algorithm (25 Da frame 
size), baseline correction by 2 runs of the multi-polygon 
algorithm, and peak search by spectra differentiation. 
The number of peaks was limited to 100 per MSP and 
all peaks were normalized to the most intense peak with 
an intensity of 1.0. The minimum frequency of occur- 
rence within the 24 single spectra was set to 50% for 
every mass. Peak lists of MSP were exported for further 
evaluation. 
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